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We derive transport equations for fermionic systems with a space-time dependent 
mass matrix in flavor space allowing for complex elements leading to CP violation 
required for electroweak baryogenesis. By constructing appropriate projectors in fla- 
vor space of the relevant tree level Kadanoff-Baym equations, we split the constraint 
■ equations into "diagonal" and "transversal" parts in flavor space, and show that 
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they decouple. While the diagonal densities exhibit standard dispersion relations at 



Q_i' leading order in gradients, the transverse densities exhibit a novel on-shell structure. 

Next, the kinetic equations are considered to second order in gradients and the CP- 
violating source terms are isolated. This requires a thorough discussion of a flavor 

x. 

independent definition of charge-parity symmetry operation. To make a link with 
baryogenesis in the supersymmetric extension of the Standard Model, we construct 
the Green functions for the leading order kinetic operator and solve the kinetic equa- 
tions for two mixing fermions (charginos). We take account of flavor blind damping, 
and consider the cases of inefficient and moderate diffusion. The resulting densities 
are the CP-violating chargino currents that source baryogenesis. 



PACS numbers: 98.80.Cq, ll.30.Er, 11.30.Fs 



* T.Konstandin@thphys.uni-heidelberg.de T.Prokopec@phys.uu.nl; M.G.Schmidt@thphys.uni-heidelberg.de 



2 



I. INTRODUCTION 

The problem of fermion mixing in kinetic theory is important both for electroweak scale 
baryogenesis, as well as for neutrinos. In the case of baryogenesis one deals with the dy- 
namics of chiral fermions at the electroweak phase transition, with a mass matrix generated 
by a space-time dependent Higgs condensate. For studies of baryogenesis in the Minimal 
Standard Model the relevant mixing occurs between the standard model quarks QJ^LFot 
baryogenesis in two Higgs doublet models the relevant mixing is in the Higgs sector 




In supersymmetric extensions o: 



neutralino sectors 



insions of the Standarc 

y,y,H 11,0,0, 3 



Model the mixing is in the chargino and 

3,0, Oi. 



We consider the dynamics of chiral fermions interacting with a plasma according to the 
Lagrangian, 

C f = itp fill) - i> L mil) R - 4> R m*4> L + £ int . (1) 

where m = + im a is the fermion mass matrix and C- m t specifies interactions of fermions 
with the rest of the plasma. Since we are interested in a near equilibrium dynamics of 
fermions at the phase transition where they acquire the mass through a Higgs mechanism, 
we shall assume that the mass is space-time dependent. Moreover, we shall assume that a 
'thick wall' approximation applies, in the sense that the typical momenta of particles are 
much larger than the rate of change of the background, d^m <C k^m. Since we are mostly 
interested in electroweak scale baryogenesis, we shall assume that £ f violates CP symmetry 
either through complex elements of m^, or through interactions C mt (complex Yukawas, 
etc.). 

Due to space-time dependence of the mass matrix, CP violation, which is crucial for 
baryogenesis, is present already in the case of two fermion mixing. 

Basis of our discussion as in Refs. }l^.ll7| are the Kadanoff-Baym equations for fermionic 
Wightman functions. They are Wigner transformed such that one can eventually detect 
semiclassical features in the resulting transport equations. In the case of several mixing 
flavors, semiclassical (quasiparticle) dynamics does not necessarily lead to an accurate de- 
scription of the kinetics. In this work we investigate the kinetics of mixing fermions in 
gradient approximation, but without resorting to a semiclassical approximation. 

In Refs. jld Ir^ such quasiclassical behavior was found at order h in gradient expansion. 
The semiclassical behavior was argued on the basis of a derivative expansion and a separation 
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of the basic equations into constraint and transport equations. The latter turned out to be 
naturally consistent and allowed reduction to transport equations for the spin dependent on 
shell distribution functions in the position dependent mass eigenbasis. Exept in the case of 
a near mass degeneracy, the dynamics of nondiagonal elements of the Wightman function 
could be shown not to influence the dynamics of the physically interesting diagonal ones 
to order H of the derivative expansion. In this work we relax this limitation combining the 
basic equations differently and we are able to discuss the nondiagonal Wightman functions 
and their highly nontrivial constraint equations, which in general are not algebraic. The 
resulting formalism does not depend on a particular choice of basis, resolving thus the 
principal limitation of the formalism presented in Refs. [lfi. 17|. 

A similar attempt within the same Schwinger-Keldysh formalism has been made in 



12j|. The principal limitation of that work is that the projection onto the diagonal densities is 
made before the relevant charges are allowed to propagate, making thus transport of mixing 
fermions unfeasible. Moreover, the CP-violating currents are inserted as sources into the 
transport diffusion equations in a phenomenological manner, which does not come out of 
the Schwinger-Keldysh formalism. 

This work bypasses both of these limitations, albeit in a somewhat simplistic disguise. 
Namely, we assume damping to be flavor blind and momentum independent, having as a 
consequence the following two limitations. First, we cannot model different damping rates of 
diagonal and off-diagonal densities, which may be of crucial importance for proper tracking 
of flavor decoherence. Second, taking a momentum independent damping may give a naive 
picture of transport (diffusion) of CP-violating currents. These limitations are of technical 
rather than fundamental nature however, and can be overcome by a more fundamental 
treatment of collisions, which is the subject of a forthcoming publication. 



II. TRANSFORMATION TO THE CHIRAL SYSTEM 

In order to study the near equilibrium dynamics of fermions in the presence of flavor 
mixing, it is convenient to start with the Kadanoff-Baym equations for mixing fermions, 
which can be derived from the effective action in the Schwinger-Keldysh formalism jlfij] . 
When written in Wigner space, the Kadanoff-Baym equations for the Wightman functions 
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iS ,<,> become (with flavor and spinor indices suppressed): 

(* + % - p -m h e-i° d " - i^m a e-^ d ^S<> - e-^{Z h }{S<>>} - e^E^}^} = , (2) 

where S h = S* - (S> + S<)/2, E h = E* - (£> + E<)/2 (5* and E* denote the time ordered 
(chronological, Feynman) Green function and the corresponding self-energy) and the collision 
term reads 

Cip = ie-*({E>}{£<}-{E<}{£>}), (3) 

with o{a}{6} = (l/2)(d x a) ■ d^b — (l/2)(^a) • d x b. These equations are formally exact, 
provided the self-energies E <>> , E a = (z/2)(E > — E < ) = and E^ are calculated exactly. 
(Note that E^ is also given in terms of E <,;> through a spectral relation ^(|). Usually one 
is lead to a reasonable approximation scheme for the self-energies. The one often used is 
based on a truncation of self-energies at a certain loop order (for an example of a one-loop 
calculation of the self-energies see Ref. Q). Equations ((2HHJ) fully specify the dynamics 
of fermions, since Sh can be determined from S <,y through a spectral relation. In some 
situations a more convenient system of equations may be the equation for S K (or F = 
(1/2)(S < + S > )) together with the equation for the spectral function A = (i/2)(S > — S' < ), 
which is collisionless, 

(k+^p -m h e-^-^ - i^ mae -fd-d k ^ A _ e -io {Eh}{A} _ e -io{T f }{S h } = , (4) 

where rrih = (1/2) (m + rrv) and m a = (l/2z)(m — rw) denote the hermitean and antiher- 
mitean parts of the mass matrix m, respectively. Now we assume that C- m t is governed by a 
set of weak couplings, and that we are interested in nearly equilibrium dynamics of the modes 
whose momenta and energies are not very small (when compared with the temperature of the 
system), such that the perturbative approach is justified. In this approach one first solves 
Eqs. (QHH) for the thermal tree- level propagators (usually in a gradient expansion), and then 
uses these propagators to study the near equilibrium dynamics by recasting the equations in 
a linear response approximation, with suitably truncated self-energies. This approach was 



pioneered in Refs. 
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12 1 . In this work we focus mostly on finding an approximate solution 
to the tree level dynamics of mixing fermions in a gradient expansion. We shall not make 
any assumption concerning the eigenvalues of the mass matrix, such that our approach is 
valid also when there are nearly degenerate mass eigenvalues, which is not the case with the 



approach advocated in Ref. Q|, which applies when h 2 k ■ d <C mf — m| ( Wi 7^ j), where mf 
are the eigenvalues of the mass matrices squared mm) and rrvm. 

For simplicity we assume planar symmetry, such that in the wall frame m = m(z), where 
z denotes the direction along which the wall propagates. This assumption is justified when 
the bubbles become sufficiently large 3]. Working in this frame it is not hard to show that 
the tree level Dirac kinetic operator V defined by the equation (c/. Eq. ©) 

ViS^ = (fc+^@- m h e? dzdk * - i7 B m a e* a * % «)iS ,< ' > = , (5) 

commutes with the spin operator, 

[V,S,\ = 0, S z = i(/A; -7 1 ^-7 2 ^)7 3 7 5 , h = sign(fc )(A; 2 -^-^ ) i/2 (g) 

k 

provided the coordinate dependences of the Wightman functions are of the form iS K,> = 
z5 ,<,> (k, t — (k x x + k y y) j k , z) . In the rest frame of particles S z measures spin in ^-direction, 

rest frame 

such that S z > 777. Having found a conserved quantity, we can write the solution 

of (0) in a block-diagonal form in spinor space (diagonal in spin) 

iS<>> = ^<'> , iS<>> = -P. [S^ g ^> _ S7 3^<,> + ^<,> _ i7 6^<,>] j (7) 

s=±l 

where P s = |(1 + sS z ) is the spin projector, P S P S > = S SS 'P S (s,s' = —1,1), and g^, g*, 
g\ and g\ denote vector, scalar, pseudo-scalar, and pseudo-vector densities of spin s on 
eight-dimensional phase space {k,x}, respectively. 

Upon multiplying Eq. (j5J) by {P s t, P s j , —P s isjs, —P s js} and tracing over spinor one 
finds Q 

f 2i k - kodt + kvVll )g s - (2isk z + sd z ) g\ - 2im h e&*>gl - 2im a e^g s 2 = (8) 
v k / 

f 2 ik - kodt + k W V \\ \ g s _ ( 2skz _ lsdz) g s _ 2im h e^g s + 2m a e&**g' 3 = (9) 
v k J 

f 2i k - fe od* + fc,|-V|| X . + ^ _ ^ g s _ 2mhe id z aZ g s _ 2im a e&£>g' = (10) 
v k ' 

(2i~k - k ° dt + fc " ' Vll )g s 3 - (2isk z + sd z ) g s + 2m h e^g s 2 - 2m a e&**g' 1 = 0, (11) 

where k» ■ Vy = k x d x + k y d y , and we have dropped the superscripts <, > of g s a (a = 0, 1, 2, 3). 

This basis is useful in the one flavor case (as well as in the mixing case with well separated 
mass eigenvalues), since at order h in gradient expansion the vector density g obeys an 
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algebraic constraint equation, from which one obtains the dispersion relation with a spin 
dependent CP-violating shift appearing at order h. An important implication of this result 
is that, at order h, the quasiparticle picture of the plasma is preserved j3, Q. When 
inserted into the kinetic equation for g^, and upon integration over the positive and negative 
frequencies, one arrives at the Boltzmann-like kinetic equation for the distribution function 
for particles and antiparticles, respectively, which at second order in gradients (first order 
in h) exhibits a spin-dependent CP-violating force. 

In the case of several flavors however, the basis g s a leads to mixing between different <^'s 
already at the classical (leading order) level. Moreover, g^s do not transform in a definite 
manner under flavor rotations [16]. A more appropriate basis to describe fermion mixing is 
the chiral basis 

g s R = g s o + g s 3, 91 = 90-9%, ( 12 ) 

and the following densities, 

g s N = 9t + W2 9 s n = 9 s i - i92 ■ (13) 
These densities do transform in a definite way under mass diagonalization (flavor rotation), 

m — > rrid = UmV^ , m) — > m d = Vrn^U^ = , 
mm) — > m 2 d = Umm^U^ , rrvm — > m 2 d = Vrn^rnV^ , (14) 

where the unitary transformation matrices U and V are chosen such that = mJ d are 
diagonal mass matrices with real eigenvalues 



9n - 9m = Ug s N V\ gtf -> g s N J = Vgtfrf. (15) 



From Eqs. (J8 HTT|) we easily find 



k 
■k\ 



= Ug s L U\ g s R - g s Rd 


= Vg s L V\ 




= Ug s N V\ gtf -> g s J 


= Vg s N ] U^. 




-)g s R - s(2ik z + d z )g s R - 


- 2irn f Eg s N 


= 


T 

+ s (2ik z + d z ) g s L - 


- 2imEg s N ] 


= 


^jg s N + s {2ik z + d z ) g s N 


— 2imEg s R 


= 


)g s N ] -s(2ik z + d z )g s N ] 


- 2im ] Eg s L 


= 



2 fk - kodt + kvVll )g R - S [2ik z + d z ) g s R - 2im)Eg s N = (16) 



2ifc - 11 !L )g s L + s (2ik z + d z ) g s L - 2imEg 3 ^ = (17) 

2i ~k - k ° 9t + r fc " i V " )g' N + s (2ik z + d z ) g s N - 2imEg s R = ( I S j 



2l ~ ko - kodt + k W V \\ \ g M _ s ( 2ikz + d z ) g^ - 2im?Eg s L = , (19) 
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where we introduced the following notation 

E = exp(~0,d fc ,) , £t = exp ( _ ~d kz d z ) . (20) 

Definite transformation properties of these equations are apparent. Indeed, Eqs. (fTHHlTjj) 
transform just as the densities g s Ll g s R , g s N and g s N * in Eq. (fTB]). 

From the antihermitean parts of Eqs. (|ToT[T7|) we get the corresponding constraint equa- 
tions for g s R and g s L , while the constraint equation for g s N is obtained simply by taking a 
hermitean conjugate of (fl"9|) and subtracting the result from |T8|) . 

(2k - 2sk z )g s R - m^Eg s N - g s N ^E^m = (21) 

(2k + 2sk z )g s L - mEg s N ] - g s N EW = (22) 

(2~k - isd z )g s N - mEg R - g s L &m = (23) 

(2k + isd z )g s N ] - m)Eg s L - g s R &rrJ = . (24) 



Note that the constraint equation for g s N ^ is simply a hermitean conjugate of 

Analogously, the kinetic equations for g s R and g s L are obtained from the hermitean parts of 
Eqs. fll6HfT7j) . while the kinetic equation for g s N is obtained by adding a hermitean conjugate 
of JEIJ) to ((TED, 



k d t + k\i ■ V| 
ko 

k d t + fcii • Vii 



ko 

k d t + fen ■ Vi 



^ - sd z g s R - inJEg s N + igtf&m = (25) 
^ + sd z g s L - imS^ + ig s N EW = (26) 



fc 

k d t + feii • Vi 



+ 2zs/c^ - imigj + ig\pm = (27) 

■gtf - 2isk z g^ - irrJEgl + ig s R &rrJ = (28) 
k 

As above, the equation for g^ is a hermitean conjugate of ()27j) . The collision terms and 
the self-energies of the constraint (|2*lTl2^j) and kinetic equations (|23H2^j) can be easily re- 
constructed from Eqs. (j2H3) by using the methods developed in Refs. r], Il2|, and shall 
be addressed elsewhere. The kinetic and constraint equations represent an exact tree-level 
description of fermionic dynamics in the presence of bubble walls with planar symmetry. 

We will now show how to solve equations (|2~lTl2*4"j) and (J2SH2HJ) in a gradient expansion, but 
without performing flavor rotations before we decouple partially the equations, in contrast 
to what was done in Refs. jlfl 17 1. 
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III. CONSTRAINT EQUATIONS 

To get an idea about the classical quasiparticle limit of our solutions, we first consider the 
constraint equations to lowest (classical) order, which are obtained from (|2~lTl24)) by taking 
the limit E ->• 1 and -»• 1, 

2(k - sk z )g R - w)g N - g ] N m = (29) 
2(k + sk z )g L -mg J N - g N m ] = (30) 
2k g N - isd z g N - mg R - g L m = , (31) 

where for notational simplicity here and in the subsequent text we drop the superscript spin 
index s. In order to solve these equations to lowest order, it is convenient to make use of 
the self-consistency of the system of equations (|2*TH2^j) and (J2SH2EJ) (c/. Ref. |3]), and use 
the solution of the kinetic equations (|27M28|) to lowest order, and working in the stationary 
limit, in which g R = g R (k^,z), gi = giik^.z) and = g^ik^jz), and where d t and Vy 
derivatives vanish, 

9n = 7rr( m 9R- 91™) (32) 



2sk 
1 

2sk_ 

When these equations are inserted in (J2DHSDI), one gets 

(k - sk z - —-{ m ^m, -}^g R + l^ m] 9Lm = (34) 
(ko + sk z + -^-{mm\ -}^g L - ^r m 9R m ^ = 0, (35) 

where {a, b} = ab + ba denotes the anticommutator. These equations can be decoupled by 
multiplying (|3*3jl by m from the right and by m* from the left, and insering the solution 
into (|3*%jl (and performing an analogous procedure for the other equation). The result is 

lr ♦ , 1 



W = 77V (gR m ^ - m] 9L ) • (33) 



(k 2 - -{mW,-} - m ] m,[m ] rri,-} ^ g R = (36) 

^ 2 -^{mmV}-^ mm\[mm\-} ^jg L = 0, (37) 



where k 2 = k\ — k 2 z = k\ — k 2 , and we made use of {a, {a, /}} — 4a fa = [a, [a, /]] . These 
constraint equations are easily solved by transforming to the diagonal basis, that is by 
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applying on the first (second) equation V (U) from the left and {U^) from the right, we 
find that the mass shells of and Ql are identical, 

{k 2 + m)) - ^ K 2 - m )f) (9R/L d h = (38) 



where mf (i = 1, .., N) are the (diagonal) entries of the matrix m 2 d = diag(m 2 , m|, .., m 
VrrvmV' = UmrwW. The solution is given by the spectral form, 



{VgfV% = 27r(fc + S A; 2 )sign(A;o)5(fc 2 -^(m2 + m J 2 )-^K 2 -m J 2 ) 2 )^.(fc,x) (39) 
(Ug s <U% = 2ix(k Q - sk z )sign(k )5 (k 2 - ±(m? + m 2 ) - ^(m 2 - m 2 ) 2 )r^(A;, x) , (40) 



where nf - = nfj(k,x) denote distribution functions. (The indices s and < are here restored 
and k = sign(k )(kl — k 2 x — k 2 ) 1 ^ 2 .) It now immediately follows that the dispersion relations 
for the densities nf • are given by 

±k = = (k 2 + \{m 2 + m 2 ) + ^(m 2 - m 2 ) 2 )^ (41) 

From (J3*9TI4*U|) one infers that, while the diagonal densities are projected on the standard 
classical shells, uju = (k 2 + m 2 ) 1 / 2 , the shells of the off-diagonal densities nfj (i ^ j) are 
given by Uij in (|4*T|l . and are in principle different for each choice of i, j. 

A particularly simple case is when there are only two fermionic flavors and only one 
off-diagonal shell. In this case gi and gn can be decomposed into diagonal and transverse 
densities as follows, 

9r = 9r+9 T r, 9l = 9E + 9l, ( 42 ) 

where 

9r = Pr9r, 9E = PF9li 9l = Pl9R, 9 T l =P T l 9l. (43) 



The projection operators are defined as 

P R = [™ + ™> •]] > P R = 1 - P R ( 44 ) 

P l = l^yi [ mm ^ l mm ^ •]] > P l = 1 - P l (45) 

where (m 2 ) 2 = [tr(m^m)/2] 2 — det(m^m) = [tr(mm^)/2] 2 — det(mm'). This notation has 
its origin in rewriting mm) = Yl 3 a=o( mm ^) a °' a m terms of the Pauli algebra, a a = (t,a l ), 
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[a 1 , o"- 7 ] = 2ie^ l a l , such that m 2 = X^Li( mm ^) 1 = Si=i( m ^ m )% where we have used the fact 
that m 2 is an invariant. In the frame in which m)m (mm) ) is purely diagonal, g R (g®) 
are diagonal, and g\ (g^) are off-diagonal and thus transverse, which explains the notation 
in (IIHH. 

Using the projectors Pr/ l and Pr/ l , the constraints for the diagonal and transversal 
parts decouple and read 

(k 2 -mm ] )g1 = (46) 
(^ 2 -M°-lf)^ = (47) 

In deriving these equations we used, [mm) , g s L D ] = 0, {mm) , g^} = tr{mm))gf j = 2(mm')° 'gj 
and [mm\ [mm\g r £\\ = An^gJ^. Similarly, we have 

(k 2 -m)m)g% = (48) 
{k 2 -{m^mf-^-)gl = (49) 

Since both (m)m) and mm) are hermitean, [m [ mf > = (mm')°, and the dispersion relations 
for the L- and R- chiralities are identical at the leading order in gradients, 

of = = = + {UmmW)«f' 2 (50) 
u T = ul = u; T R = (k 2 + (mrrJf + ^) ^ (51) 



An analysis of the constraint equations ()21H24|) shows that at higher order in gradients the 
diagonal and transverse shells mix in a manner which includes the derivative dk z , leading 
to nonalgebraic constraints for the Wightman functions, and thus seemingly breaking the 
quasiparticle picture of the plasma, which questions the validity of any on-shell description 
of the dynamics of CP-violating densities, which necessarily involve higher order gradients. 
The situation is more complex however, than this simple argument seems to indicate. As 
we show in the next section, in spite of this problem with the constraint equations, one can 
solve the tree-level kinetic equations to an arbitrary high order in gradients, thanks to the 
fact that, in stationary situations, the kinetic equations (J25M28)) do not involve ko, and thus 
the tree-level dynamics of the Wightman functions gn, gi and gx and the corresponding 
on-shell densities (obtained by /co-integration) are identical. (In nonstationary situations off- 
shell effects may be important however, which is indicated by the k dependences appearing 
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in the d t and Vii derivatives in Eqs. (|25H28jl .) This also means that stationary tree-level 
dynamics is completely specified by the on-shell solution of the corresponding Dirac equation, 
which is by no means true in general situations. 

The importance of the leading order analysis of the constraint equations presented here 
stems from the fact that it allows for the on-shell projection of the collision term and self- 
energies at leading order in gradients, such that it is essential for a self-consistent derivation 
of the kinetic equations for mixing fermions, provided one approximates the collision term 
at leading order in gradients. 



IV. KINETIC EQUATIONS TO LOWEST ORDER 
Using (jS2HSSI) m (EH) and (f2*U|) to lowest order, and working in the stationary limit, we 

get 

dz9R + T^r [ mtm ' 9r] = (52) 
d z 9L + [mm\ g L ] =0. (53) 

From the solutions of these equations 

gn(k z ,z) ^ exp (~2F^ dz'[rrJrn{z'),-]}g R (k z ,$) , (54) 
gi(k z ,z) ~ exp(-^-^ dz'[mm^(z'),-}^g L (k z ,0) (55) 

we see that the diagonal and off-diagonal densities (when viewed in the diagonal basis, = 
VggY\ gid = UgiU^) exhibit a qualitatively different behavior. The diagonal densities do 
not evolve, while the off-diagonals exhibit the vacuum oscillations, well known from the 
neutrino studies. Note the identical evolution of the L and R chiralities, when viewed in the 
diagonal basis. Note further that in the case of two mixing fermions, we have [rrvm, g^] = 0, 
[mrrv,g^] = 0, [rrvm, g^\ = 2i(m 2 x g^) • a and [mm^g^] = 2i(m 2 x g£) • a, such that the 
transverse densities rotate with the frequency specified by Co = m 2 /k z . 
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V. KINETIC EQUATIONS TO SECOND ORDER 

In stationary situations, one can rewrite the system of kinetic equations ()25H28p in terms 
of the chiral densities g R and and gi only, valid for all orders in gradient expansion, 

d z g R + ~(^E^(mEg R )) - l -(^(g R &rrf)&rn>) 

- % -(^E^(g L &m)) + % -(^-{m)Eg L )&m>j = (56) 
d z 9L + % -(mE^(m^Eg L )^ - % - (j-{g L &m)EW ) 

-~(mE—(g R &m^ + ^(—(mEg R )&m^ =0.(57) 

Note first that the chiral densities g R and gi couple through derivative terms only, which 
justifies the use of the chiral densities in writing the kinetic equations for mixing fermions. 
Next, equations ()56|) and (JSTj) are transformed into each other by the following replacements, 
R <-> L, m <-> rw and s «-> — s (see e.g. Eqs. ()34H35p ). defining thus the symmetry, 
which relates the dynamics of the chiral densities g s R to g s L . Furthermore, we have arrived 
at Eqs. ([56fJ57J) without using the constraint equations (|2"IH%iJ) . This procedure has the 
advantage that fco appears nowhere in Eqs. (|5fitiH7j) . implying that the kinetic equations for 
the distribution functions f Rs ± and /l s ±, defined as the (positive and negative) frequency 
integrals of g s R and g s L , have exactly the same form as (|56H57|) . resolving thus the problem of 
closure of the on-shell kinetic equations. We emphasize that the (tree-level) closure is thus 
achieved, even though the constraint equations are nonalgebraic. One consequence of the 
nonalgebraic nature of the constraint equations is a coupling between the off-diagonal and 
diagonal densities, which is nevertheless implemented in a self-consistent manner into the 
kinetic equations (j56tiBT|) (through the higher derivative terms), without ever referring to 
the on-shell structure of the system. If one had attempted to further decouple the equations 
for g R and g^, one would have found out that this could be achieved by making use of the 
constraint equations (J2TH2H) . which would reintroduce the dependences on k and s, which 
is not explicit in equations (JKEHKTj) . 

Upon expanding E and E^ in (|27]|) to second order in gradients, 

E = l+ l -d z d kz -^(d z d kz ) 2 + .. 

E ] = l-^d kz d z -^(d kz d z ) 2 + .. (58) 
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we can write the chiral kinetic equations (|56H57|) . truncated at second order in gradients as 
follows, 

k z d z g R + - [m f m, g R ] 

1 1 / 1 / 

--{(m t m) / \d kz g R } + ^r( mt mg R + g R m ] m') - ^(^ f 9Lm + m ] g L m!) 

[(m+m)", dl z g R ] + [m+V, d kz g R ] 
+^"md kz { 9 f )-d kz (f)rrJm") -i(rn)"d kz {f )m-m% z (^)m") ^ . (59) 

The kinetic equation for gi is obtained from (|59|) simply by exacting the replacements, 
9r ^ 9l and m rrv . 

To get a rough idea on what are the criteria for the applicability of the gradient expansion, 
let us first recall the relevant criterion for the one fermion case, 

W < hk» . (60) 

Since d z ~ 1/L W , this criterion was used to coin the term "thick wall regime" in baryo- 
genesis studies. The proper interpretation of f|UU|) is closely related to the validity of the 
WKB approximation in quantum mechanics, which is valid when the de Broglie wavelength 
of the excitations is small in comparison to the region over which the background (mass) 
varies. In the case of mixing fermions however, an additional complication arises from the 
evolution of the off-diagonal densities. Indeed, from the leading order solution (jB3|) we find, 
dk z (g Rt L)ij ~ —i[ Jq dz'{mf — m])/ (Zk^^gx^ij, such that the dynamics of the transverse 
densities can jeopardize the validity of the gradient expansion. Indeed, the gradient expan- 
sion applies provided formally, [J* dz'(mf — m^)/{2k 2 z )\ d z 1 (Vi,j) is satisfied, such that 
the dependence on L„ ~ ~ J Q 2 dz' roughly cancels out. This implies an additional and 
qualitatively new criterion for the validity of the gradient expansion, mf — m| k 2 z (Vi, j). 

This criterion should be taken with great caution, however. Being derived without any 
reference to decoherence of off-diagonal densities, this criterion may be in many situations 
too stringent. In particular, large differences in mass eigenvalues and low momenta k z imply 
fast oscillations of transverse densities, which are more prone to decoherence by rescatterings, 
and thus destruction, than slowly oscillating densities. We thus conclude that a complete 
analysis of applicability of the gradient expansion for mixing fermions is at the moment not 
available. 
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A. Reduction to the diagonal limit 

Let us now consider the diagonal part of the kinetic equation (|59|) and its left-handed 
counterpart, 

k z d z g R - ^{{m ] m)' D ,d kz g R } + ^{(m f m)'^ g R - g L ) 



{ {rn)'m — m)m')' , d k — — } = (61) 

16 I Ky J 



k z d z g L - ^{{mm ] )' D ,d kz g L } + ^-{(mm r )' D ,g L - g R ) 

where g R d = Vg R V^ = Ug R W and g^d = UgiU^ = Vc/lV^ are assumed to be diagonal. 
Since in the diagonal basis {V{m)m )'yt) = {ml) 1 = (U{mrriyW) , and 

(y{m)' m — m)m")V^ D = (U{raw)" — m"m))U^ D , (63) 

we see that the equations for g od = (g Rd - g L d)/2 and g 3d = (g Rd + gLd)/2 decouple, 

k z d z g od - -{m 2 d )' 'd kz g od - ^-(V(m)'m - m)m')'V ] ) D d kz g Qd = (64) 

dn9u-\{mWdj£ = 0. (65) 

Note that the leading order solution, gfj = (sk z /k )g^J solves (JfiHj) . where g^J is given 
in (J86|) . such that up to second order in gradients there is no source for the axial density in 
the diagonal approximation. 

On the other hand, the form of the vector equation ()64|) suggests that in the static case 
can be solved exactly. This is in fact not quite so, since the quantities in (p)3*j) are not a 



total derivative. Nevertheless, an approximate solution can be found for a static wall 16, 17]: 

g od ~ 27r\k \ d{k 2 - \m d \ 2 + ^-{V{w)'m - m ] m')V ] ) D ^jn (k ) , n = - - . (66) 

This is a good approximation when the mass matrix is approximately diagonal, or when the 
transverse elements of VV^' and UW are small. Eq. (|64|) reproduces the result first derived 
in Refs. 13, 



14 1, where the last term sources a CP- violating current, of crucial importance for 
electroweak scale baryogenesis studies. In sections lVlll IVllTl and UXl we make a quantitative 
comparison of the second order diagonal source in (p>f|) with the first order (transverse) 
sources, which we discuss in detail in the next section. 
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VI. CP- VIOLATING SOURCES 

Let us define CP symmetry as the following transformations of the Dirac spinors (up to 
an irrelevant phase), 

*fj cp {u) = CVtp{u){CVy = i^ T {u) , ij cp (u) = CV^(u){CV) ] = *fj T {u)i-f 2 (67) 

One finds that the kinetic equation (j2J) transforms as [16], 

(jt+ -fa -ml{x)e-&*- 8 * +i^m* a (x)e-^-^S cp< ' > (k,x) = - 7 ° 7 2 C;(-fc, x)-f-f , (68) 

where we neglected the self-energies. The Wightman functions transform as, 

S < ' > (k, x) ^ - 1 2 S > > <T (-k, x) 7 2 = S cp< ' > (k, x) (69) 

A comparison of Eqs. (J2J) and ()68|) reveals that in the wall frame, in which m^ a — m h,a{z), 
the CP transformation of the flow term is in our context equivalent to the transformations 

rrih — > m* h , m a — > —tn* , (m — > m*, m} — > m T ) , (70) 

leaving e.g. k~^ x d z and d z d kz invariant. From these rules the CP transformed equation 
can be written as follows, 

kAg% + i[m T m\g*\ (71) 

(i T *\i o cpi 1 / T 1 * cp , cp T *'\ 1 I T' cp * , T cp */\ 

--{(m m ) ,d kz g£\ + —{m mg£ + g£m m ) ~ Jj^{ m 9h m + m 9l m ) 



z 

I 



-± [(m T m*)" , dlg R ] + ^ [m T 'm*', d kz g R ] 
+ l -(m T "m*d kz { g f) - d kz { 9 f)mV) - 1 (m T "d k , (f)m? - m T d kz {f)m*") ~ 0. 

Our primary goal is to identify the CP-violating sources in the system. A naive way of doing 
that would be to subtract Eq. (fTTj) from Eq. (jo^j) . and identify the terms in the equation 
for 5g c R p = gR — g c R p which involve a CP- violating operator acting on g R = (g R + g R )/2, thus 
representing mixing of CP-odd and CP-even densities. This procedure leads to equations 
with indefinite transformation properties under flavor rotations, which is a consequence of 
the indefinite transformation properties of the newly defined densities 5g c R p and g R , making 
it difficult to disentangle the genuine CP- violating densities from the apparent, but possibly 
spurious, CP-violating densities. 
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In the following, we propose a method, which allows to extract the CP violation from 
the kinetic equations (JoTJj) and ([71)1 . We see from Eq. (jo^|) not only that CP symmetry is 
broken by the theory, but also that it depends on the basis in which the usual definition 
of CP conjugation is used. Therefore it is helpful, instead of the usual CP conjugation, to 
introduce another operation 

Qg(k,x) = CVg T {k,x) (C7>) f = g{-k,x). (72) 

For the second equality, the relation ((00)) has been used. 

Note that we do not want to undo a part of the usual CP conjugation, but just want to 
include an additional sign change of the imaginary elements in flavor space (the coefficients 
of o"2 in the case of two mixing fermions). This operator has some nice properties, that 
are absent in the case of the standard CP conjugation. For example, Q commutes with 
basis transformations, implying that if a quantity transform in a definite manner under Q 
conjugation in the mass eigenbasis, it will transform in the same way in the flavor basis. 
(This can be easily checked by noting that CP V (CV)^ = V*,CVU (CP) t = U* .) In addition, 
Q conjugation agrees with the usual CP conjugation for the remaining coefficients in flavor 
space (of 1, o"i and 03 in the two flavor case); in particular it agrees for the diagonal elements. 
Hence, in order to generate a CP-violating effect, one has to create Q-odd terms at least in 
the off-diagonal terms in the mass eigenbasis, such that CP violation becomes manifest in 
the diagonal terms of the flavor basis. 

Knowing this, we look for the transformation properties of the kinetic equation under 
the Q conjugation and observe that all terms that come from an even order of the gradient 
expansion acquire an additional minus sign relative to the odd contributions from the gradi- 
ent expansion. As a consequence, in the one flavor case one has to take second order terms 
into account in order to break Q (in this case Q is of course equivalent to CP), since there 
are no zeroth order terms. This leads to semiclassical force induced baryogenesis. 

In the multi flavor case, an expansion up to first order is sufficient as long as the zeroth 
order contributes (otherwise the Green function is everywhere diagonal in the mass eigenbasis 
and the problem reduces to the one flavor case as discussed in Q ) • 
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Let us recall the kinetic equation of the right handed density to first order 

i 1 
k z d z g R + - [m^m,g R ] - -{(m f m)', d kz g R ) 

+ — — (rn) mg R + g R m^m') — {m)' g^ra + m^iUi') ~ 0. (73) 

The Q-conjugate equation is 

k zd z g R - Jpm,^] - ]{{m^m)',d kz g^} 



+ 



— (rr^'mg^ + g®m^m'^ — (m)' g®m + m^g^m') ~ 0. (74) 



4fc 2 v w y 4fc, 

such that only the sign of the second (zeroth order commutator) term is affected. 

To solve this equation, we will first determine the lowest order solution and then expand 
around it. The best way to determine it is in the mass eigenbasis, since in this basis the 
'direction' of the mass in flavor is fixed. In the flavor basis we can mimic this property by 
adding a term that explicitly compensates for the ^-dependent basis transformation. In the 
case of a static wall profile, we can in addition include the diagonal part of the third term, 
since it belongs to the classical Boltzmann-like flow, and we know how to handle it from 
the one flavor case. Diagonal means in this context that it commutes with rrvm as it is 
indicated by our notation introduced in section 11111 Our lowest order solution fulfills 

k z d z gf + i[mtm,sg>] - \{{rrJm)>, d kz gf} D - k z [v+V,^] = , (75) 

where the first commutator term vanishes, since there is no source for the transversal parts. 
Now since V{(m^m)', dk^g^^V^ = 2(m^)'dk z g R l, the solution of (|73j) is simply, 

9 f =(l + SW°V. (76) 
v fc y 

Here g^J is the diagonal vector density in the mass-eigenbasis of the spectral form, 

g$ = 27c\~k \S(k 2 -\m d \ 2 )n, (77) 

and n = n(kok\i) is a distribution function. In thermal equilibrium, which is formally 
obtained in the limit of large damping (frequent collisions), n reduces to the Bose-Einstein 
distribution, n — > n = l/[exp(/c /7 1 ) + 1]. In this case the collision term vanishes (this is also 
obtained by imposing the Kubo-Martin-Schwinger condition on the Wightman functions). 
All influences of the changing background are then negligible and the Green function depends 
only locally on the mass. 
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The transverse part of the deviation from g R is in the next order given by 



Kd z9 T + \[nJm,gT] + W g™ - k z [v*V,g™ 



S 



(i) 

R 



s 



(1) 

R 



k z 
1 



f t ' (°) 

m' m — m'm ,g R 



~(o) 
91 



9r] 

1 T 1 

8kz 



' t V (°) 

v m'm) ,g R 



40) V 

9l) 



(78) 



where we defined g L = V^Ug L WV. Here we have introduced a damping rate T (which 
can be arbitrary small) to fulfill boundary conditions at infinity (g R — > for z — > ±oo). At 
the same time this helps to cure the infrared divergencies in the sources. For simplicity, in 
this work we assume that the damping T is flavor blind, i.e. we take it to be proportional 
to the unity matrix in flavor space. 

The (particular) solution of equation (fTHJ) in the mass eigenbasis is formally given by 

W(z,z')S^(z')dz', 



92%) 



(79) 



with the kernel 



W(z,z') = ^[9{k z k )9(z-z')-9{-k z k )9(z'-z)]e W (- 1 ^-T{ 

x ex P ( -TTT I [ m d(y)> •] d v 



z — z 



2/c 2 



(80) 



The exponential function is understood as the power series in nested commutators, and the 
source is rotated into the mass eig enbasis Sgj = V S R 1] V. 

From this we can deduce the part of g^ that breaks the Q-symmetry (g® R = g R — g R ) 



W4z,z')S^(z')dz' 
) - 9{-k z k )9{z' - z)] exp ( - ^T(z - z'\ 



with 



W (z,z') = ±[9(k z k )9(z - z 

1 



x — i) sin 



2,k z 



(81) 



(82) 



The CP- violating diagonal part of g& in the flavor basis is given by Tr(a^ Vg^^V^. 



A. Local Sources 

Since our source is already first order in gradients, we solve the integral (JHTj) by expanding 
all functions around the position z and keep only the first Taylor coefficients. This procedure 
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is justified provided r L w > 1. For the MSSM this leads to V > L" 1 ~ T c /20 Q, with the 
expansion parameter of the gradient expansion T C L W ~ 20. Since the off diagonal entries 
are coherent superpositions of particle states, F characterizes the inverse decoherence length. 
On the physical grounds we expect T to be at least as large as the thermalization rate, which 
we take to be of the order the thermalization rate for the W bosons, T ~ Tw — o. w T c 2fl ]. 
To get its detailed form and magnitude would require a quantitative analysis of the collision 
term for mixing fermions however, which is beyond the scope of this work. Here we take T 
of the order the thermalization rate and proportional to unity in flavour space. 

Assuming that T induces an efficient flavor decoupling and making a leading order ap- 
proximation of the sine function in Eq. ()82)1 . the expression contributing at leading (first) 
order in gradients acquires the following simple local form, 



T(l) 

9r 



such that the Q-breaking part reads 



9l® 



k n T 



[mtm, • ] J S 



7 (i) 
'r ' 



m^m, ■ ] 



:S 



(i) 



(83) 



(84) 



k 2 T 2 + | [mtm, [mtm, • ]] R 
Since only transverse sources lead to Q-breaking terms we rewrite the expression (|7Hj) in 
the more compact form 



S 



(i) 

R 



k 7 . 



— i ( in 

4 



4fcn 



\l f / -(0) 

m 1 m — m'm , (yfQ 



] m)\d kz gf) J 
s 

~4lJ 



t v -(o)V 



(85) 



where g% 



(0) 



V^Ug ( ° ) WV)/2, and we made use of the leading order 



constraint equation, g 3 = s(k z /k )gQ , with g^ = V^g^JV and 

gf} = 27r\k \S(k 2 - \m d \ 2 )n , 



(86) 



where in equilibrium and for a static wall the occupation number reduces to the Fermi-Dirac 
distribution, n — > n = (exp(k /T c ) + l) -1 . 

Here we can establish for the first time the important fact that in our treatment a static 
wall does not induce any CP-violating charge densities. In the local approximation (J8%j) the 
following integrals are relevant for the calculation of the CP-violating source current 

•s / d^k 7; [mtm, • ] k Q 



JR 



I [mtm, 

4 J (2vr) 4 ^r 2 +|[mtm, [mtm, -]]|g 



4ki 



v*v,§® 



+ 



\l f / -(0) 

m' m — m'm 



+ 



Wm)',^} 



(87) 
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The left-handed source current j T L is obtained simply by the substitutions, s — > —s, m <->• m) 
and g^ = V^g^Jv — > U^g^JlJ. The second term in Eq. ()85|) does not contribute to the 
current (|57j) . since the integral over the momenta vanishes for this term. 



B. Nonlocal Sources 

To evaluate (|5Tjl we could just solve the integral numerically. However this would involve 
some technical and physical shortcomings. First, the integrand is oscillating with a frequency 
uj ~ A/k z , which makes numerical evaluation hard. Second, since we have parametrized the 
collision terms in the kinetic equation by just one parameter, our solution does not show 
the expected behavior in certain regions of parameter space. E.g. we expect that collisions 
help to isotropize the deviation from equilibrium, while the solution to equation (|81|) has a 
strong k z dependence, but almost no k\\ dependence (fen denotes the momentum parallel to 
the wall). Another feature which may play an important role is diffusion, by which particles 
get transported typically to distances 

2D 

idm ~ i 1 2 i Arn\i/2 ( 88 ) 

in front of the wall, where D denotes the diffusion constant, v w the wall velocity, and T the 
damping. In systems with small T and/or large D, such that v w ^> ATD is satisfied, the 
diffusion tail may be large, £dis — D/v w . Since for charginos of the MSSM, the diffusion 
constant and the wall velocity are rather small, D ~ 10/T, v w < 0.1 and the damping quite 
large, T ~ a w T 20], ATD is amply fulfilled, and we can estimate the diffusion 'tail' to 

be fdifr — (-D/T) 1 / 2 ~ 15/T. Since diffusion is symmetric and it extends to distances of the 
order the wall thickness, we expect it to be captured reasonably well by our simple model 
of damping. 

To cure the shortcomings related to the local approximation, we shall solve the kinetic 
equation (|7S|) by a fluid Ansatz in the mass eigenbasis 

JV 

9sdij = 2it^T- a ^ aij (k z -v w k ) a [d ko no{^(k -v w k z ^ (89) 

a=0 

where are given by the lowest order on-shell conditions (jSJ), ufj = oufj — fcjj, and n (x) = 
l/[exp(x/T c ) + 1] is the Fermi-Dirac distribution function. Note that the fluid Ansatz (fKHj) 
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captures the chiral nature of the first order solution, which is, for example, expressed by the 
leading order constraint relation, <?}£L = (l + s& z /&o)#odL- 

If one now takes the first N momenta of the kinetic equation (|78)1. defined as 
§ ko>0 [d A k / '(2vr) 4 ](/c / 'ko)(k z /T c ) 1 (I = 0, ..,N), one gets a matrix equation of the form (here 
and below we suppress the i,j indices): 

Ad z fi+^B[m 2 d ,fi] +TCfi = D. (90) 

A, B, C are matrices and /x and D vectors in the a, b space (a, b G {O.JV}), with 

A ab = T- a - b / |fc |(fc 2 -^o) a ^+ 1 [9 fco n ( 7 (A;o-^A; z ))]5(fc 2 -^ 2 ) (91) 

7fc >o \ Zl1 ) 

74 Z 



B«6 = T~ a ~ b ! -^~M{k z - v w k ) a k b z [d ko n (-f(h - v w k z ))]5(k 2 - u 2 ) 

Jk o >0 \ llx ) 

C ab = T c 7 a - b f 4^\k \(k z - v w k ) a k k b z [d ko n {l(ko ~ v w k z ))]S(k 2 - to 2 



k Q >0 



x [4k 2 z [VV v ,n ] + [V{m v m - m ] m')V\n ] T + {V (m< m)'V\ n } 

The summation in Eq. (|9T?|) runs over b, while a, i and j are held fixed. 

The eigenvalues tj of the matrix A -1 ((z/2)v3[m^, •] + TC) determine the damping and 
oscillatory behavior of the solution. Due to the form of the source, the first few momenta 
dominate the solution. If the source has a compact support, we can deduce that outside this 
compact region, \i is a superposition of damped harmonic oscillations with the frequencies 
£y(7i) and damping rates 3?(7i). The amplitude of these oscillations is then suppressed by 
|7i| _1 , such that fast oscillating modes give smaller contributions to the current. 



VII. DOMINATION BY DIAGONAL PARTS 



For special choices of the mass matrix, or in the limit where the oscillations suppress the 
off-diagonal contribution, the problem can be reduced to the diagonal case and the first order 
contributions to the CP violation, produced by the oscillations of the off-diagonal terms, are 
negligibly small. When viewed in the mass eigenbasis, the problem then reduces to the 
diagonal case, such that the first CP-violating contributions come from the second order 
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semiclassical force in the kinetic equation. This approach was originally pursued in jltl Il7|. 
and we summarize its main results in section IV Al 

We have seen that the first order terms are for large damping suppressed as T~ 2 . In this 
section we pose the question how are the second order terms suppressed in this limit. 

Since T is large we integrate the Taylor expansion of the source using the Green function 
method and notice, that the first coefficient gives no contribution (since it is odd in k z ) and 
the second term gives 

(2)D isk z ( ui i „ „ -(0)V D /no\ 

9o = su2r2l l m ~ m ,dk * 9 ° i ■ ^ > 

oKqL Kq 

The term (J92|) is suppressed by T 2 as the contributions in the first order terms, but in 
addition by two more orders in the gradient expansion. Therefore we can not infer that 
the second order terms dominate for large damping. Rather the region in parameter space 
where A is large leads to dominance of the diagonal terms. However, the second order terms 
can yield CP violation in the trace of the Green function, while the first order terms are 
always traceless. Therefore the second order terms could be more important for generating 
a baryon asymmetry, depending on which contribution is more efficiently transformed into 
the BAU. 

VIII. LOCAL CONTRIBUTIONS TO THE CURRENTS IN THE MSSM 

In this section we give the explicit expressions for the CP-violating currents in the MSSM 
in the local approximation. Since in the MSSM we expect the damping T to be less than 
T c /20, we are in the regime, where transport is important, such that the main intention is to 
make our approach comparable with former publications, in which local sources for diffusion 



equations have been derived 



n 



12 



The most important contribution to the BAU in the MSSM is determined by the mass 
matrix of the chargino-higgsino sector with complex M2, fi c and real H^H^- 

m =("" HS ) (93, 

\gW ft / 

The procedure how to diagonalize m is outlined in Appendix [X] 

Using this parametrization we can evaluate the CP- violating chiral source current j^r (JHTj) 
and j T j®. Since this sources are traceless, the relevant quantities are Tr(a 3 j^) and 
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Tr (cr 3 j^), where a 3 = diag(l, — 1) in flavour space. Under the trace Eq. (J87)l can be 
reduced to the form 

is f d 4 k A kr, 



8 Jk, 



8j k0>0 (27r)3*gr«+jA!»|jfc b | 
x | -Tv[a 3 5{k 2 -\m d \ 2 )n]Ti[(Va 3 V ] ) T {Ak 2 z (VV V ) T + {V{m ] 'm-m ] m')V ] ) T }] 

+ Tr[6(k 2 -\mX)n]Tr[(Va 3 V^ T (a 3 V(m^m)'V^ T ]^. (94) 

The traces can be easily evaluated by making use of Appendix A, 

Tr[(Va 3 VY(VV^) T ] = 4*-^(M 2 /i c ) (u[u 2 - Ul u' 2 ) 

Tr[(Va 3 V^ T (V(m p m-m^m')V^ T ] = -4i^3(M 2 // c ) (uiU 2 )' 

Ti[(Va 3 V ] ) T {a 3 V{m)m)'V ] ) T ] = -U^H(M 2 H C ) (u[u 2 - Uiu' 2 ) , (95) 

where we used A = A. The form of the chiral source can be then written as (we reinsert the 
spin superscript) 

Tr(^f ) = s^^^(u lU2 y V t 0) - 3^(„' lUa - Ul u> 2 ) ( V f 0) + 4 V f 2) ) (96) 

c c c 

rp / 3-s®\ $(M 2 /i c ) A , x/ 3 . ^(M 2 fJ, c ) ( , , u o 3 x 

1 c 1 c 1 c 

where A = \M 2 \ 2 — |/i c | 2 + {u\ — uf), A = |M 2 | 2 — |/i c | 2 — (w 2 — w 2 ), and we defined the 
integrals, 

" W1/2 = i>„ WS'" *8p + (a/2)» 4( * - m " 2) 

1 . T 2 

»7(n) = 9 (^(n)l + ^(n)2j , = ^ fan)! ~ V(n)2) , A = Vtl\ - m\ . (98) 



The functions 77? % and r/ 3 ^ are dimensionless and depend only weakly on A in the region 
where \kn\T > A, but generate a behavior oc T^/A 2 in the limit T — > 0. In our former 

n 

publication [16] we neglected off- diagonals, so we were working in the limit of large A. 



Several comments are now in order. Equations (j96H97j) represent the new CP-violating 
sources which were calculated by solving iteratively the quantum kinetic equations for two 
mixing charginos of the MSSM. These sources are absent in the single fermion case, in which 
case the diagonal semiclassical force source dominates. Both of the chiral source currents 
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are proportional to spin, and hence when summed over spin the sources vanish, 

£^=°> £j?=0- (99) 

s s 

Nonvanishing sources are obtained only when a weighted sum over spin is performed (for a 
related discussion of the semiclassical force source see Refs. |3,[l2|), sj s ® and ^2 3 sj^. 
To get the source currents for a moving wall, we assume that in the plasma frame the current 
transforms as a Lorentz vector, which is reasonable provided diffusion is inefficient. In order 
to facilitate a comparison with the existing work on electroweak baryogenesis sources, it is 
instructive to calculate the vector and axial source currents. From Eqs. (jHEUHZI) we then 
easily get for the currents in the plasma frame and for a moving wall 

Tr^t) = Ito'E j (,-g + = 2 8( y W ~ !g=£ [B^Mv (100) 



2 \ J Hfl 1 •'Lip,) rp2 rp2 

s c c 

■vr^jf) s if e | at - m - 2 V "U 



7^2 7^2 

s c c 

- 2 ^f c \ u 2 d tl ui - Ul d,u 2 ) (rf( Q) + 4^ 2) ) ,(101) 

c 

such that the sources in the local approximation neatly split into the plus and minus con- 
tributions, oc d IM (uiU 2 ) and oc u 2 d il u 1 — u 1 9 A1 w 2 , respectively. The axial current is sourced by 
the plus contribution only (just like in the case of the second order semiclassical force), while 
the vector current is sourced by both plus and minus contributions. In the non-local case, 
both plus and minus terms contribute to j® and jf n . These results have a similar structure 

1 I— I l—l J M J Ofl 

to the sources found in Refs. jljj and jl2(|. The differ however, when a detailed quantitative 
comparison is made. 

Finally, we quote the second order diagonal source calculated in the local approxima- 
tion (jHS): 



Tr(ljg) c 2^^d»(u»u 2 + u lU fy? 0) (102) 



with the definitions 



C(n)l/2 = T c ^ j2^JT kz fcgp <K« - m V2) , C(n) = ^ (C(n)l ~ C(n) 2 ) • (103) 

Note that, in contrast to r/, (is not suppressed for large A, such that the second order terms 
dominate in the local regime for large values of A. 
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FIG. 1: Bubble wall: The higgs vev profile (|104|) and (3 (J105JI as a function of z. 



IX. NUMERICAL RESULTS OF TRANSPORT IN THE MSSM 



In this section we will present numerical results of the fluid ansatz. The Higgs vevs and 
the (3 angle are parametrized by Hi(z) = H(z) sin(/3(z)), H 2 (z) = H(z) cos(/3(z)) and 



H(Z): 

/3(z) = p 



Wl-tanh( a (l-£ 



^ A(3 ^1 + tanh ( a ( 1 



2z 



(104) 
(105) 



If not stated differently, the parameters used in the plots are T c = 95 GeV, v(T) = 175 GeV, 
a = §, tan(/3 00 ) = 10, T C L W = 20, F = a w T c , M 2 = 200 GeV, /i c = 250 GeV, and the 
complex phase is chosen maximally 3(M2/i c ) = \M2fi c \- The value of A/5 = 0.0108 is 
deduced from ^| by using the value = 200 GeV. In figure Qwe show how the wall ()104l - 
I105|) looks for our choice of parameters. 

The wall velocity is taken to be v w = 0.05 and the plots are evaluated for the first six 
momenta. The currents j°, j® and the second order term are evaluated in the plasma 
frame, thus the expressions are linear in v w . 

First we will display the dependence on the number of momenta that are used in our 
fluid Ansatz. Fig. El shows that the convergence is already very good for N = 6, such that 
it is sufficient to use only the first six momenta. Since the higher momenta have a slightly 
smaller exponential suppression factor but much smaller sources, these contributions become 
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FIG. 2: Dependence of the current for a typical source on the number of momenta N used. 

important in the region far away from the source, where the oscillations take place. Mostly 
the phases of the oscillations are influenced, such that quantitative statements barely change 
for N >6. 

The plots Fig. Elto Fig. Elshow all three currents for the values \i c = {205, 220, 250, 450} 
GeV and M2 = 200 GeV. For small A, respectively M2 — \i c , the solution is oscillating 
and has rather large amplitudes. These oscillations are, as expected, suppressed for larger 
values of A and a local contribution remains. For A ~ 20 T c 2 ~ T^L W the second order 
contribution, which shows a weaker dependence on A, begins to dominate. When A is large, 
the first order currents are suppressed due to efficient decoherence. In the BAU, the second 
order terms start to dominate earlier since, for the first order terms, the oscillations and 
inefficient transport prevent in part an efficient source conversion to baryon asymmetry, 
while the second order terms are transported more efficiently and without oscillations, such 
that that give a truly non-local contribution. 

The term resulting from the combination Uxd^Ui — u^d^Ux is suppressed due to the fact 
that A/3 is small. In addition, the terms that include rf are smaller than the terms including 
rf. 

The axial vector current j®, that is normally the largest contribution to the source, is 
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FIG. 3: The CP-violating currents of first order j°, j® and of second order j® . M 2 = 200 GeV. 
fi c = 205 GeV, A/T c 2 G [0.22,4.0]. 



CO 



6e- 


■06 r 


5e- 


■06 - 


4e- 


■06 - 


3e- 


■06 - 


2e- 


■06 - 


1e- 


■06 - 




- 


1e- 


■06 - 


2e- 


■06 - 




-0.4-0.3 -0.2-0.1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

z in GeV" 1 

FIG. 4: The CP-violating currents of first order j°, j® and of second order . M2 = 200 GeV. 
H C = 220 GeV, A/T c 2 G [0.93,4.2]. 
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FIG. 6: The CP-violating currents of first order j°, j® and of second order j®^ . M2 = 200 GeV. 
Mc = 450 GeV, A/T c 2 G [18.0, 19.2]. 



29 



suppressed for small A due to the factor |M 2 | 2 — |/i c | 2 in Eq. (jlUU|) and in this region the 
vector current j° becomes the most important one. 

X. DISCUSSION AND OUTLOOK 

In this article we have presented a method to solve Schwinger-Dyson equations for CP- 
violating densities for mixing fermions in a space time dependent background. The transition 
to the chiral basis was important for the partial decoupling of the different coefficients in 
spinor space. The terms that can appear in the kinetic equations written in chiral basis all 
have the same transformation properties under flavor basis transformations. As an interme- 
diate result, we have obtained the formally exact equations ()57|). in which only two of the 
16 coefficients in spinor space remain coupled to each other. 

Next, we have found that the off-diagonal densities exhibit oscillations at lowest order in 
gradient expansion. Even though they vanish without space-time dependent background, 
these terms ought to be treated by taking into account oscillations as soon as they are 
sourced. 

In section IVII we advance a novel definition of CP violation in kinetic equations with 
mixing flavors. According to our definition, when kinetic equations with mixing fermions 
are truncated at first order in gradients, only the inclusion of flavor oscillations (formally 
expressed through the commutator term) gives raise to nonvanishing CP violation in the 
diagonal entries of the flavor basis. The same is of course true for the chargino and neutralino 



his is the main difference 
lllQ- Without taking 



sectors in supersymmetric extensions of the Standard Model, 
between our approach and the approach advocated in literature 
the flavor oscillations into account, the CP- violating densities would stay in the off-diagonal 
entries even after rotation into the flavor basis. 

Our approach to second order diagonal sources (semiclassical force) differs from the treat- 



ment advocated in Refs. 



10 



16 



17| . In order to arrive at a local analytic estimate in flavor 
basis, we have considered the limit of large damping, TL ^> 1. The CP- violating axial 
current (jlU2|) is in this case proportional to the trace in flavor space. The source from semi- 



classical force in 



171 ] is calculated in the mass eigenbasis in the limit when TL <C 1, 
and it was found to be proportional to the difference of flavor axial densities, Tr(a 3 j^ d ). 
Moreover, in our numerical treatment we calculate the source by taking account of both 
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flavor mixing and transport, while the same source has been treated in the literature in the 
diagonal approximation in the mass eigenbasis. 

Apart from the plus contribution, ex uxdyUi + u^d^ux = d fl (uiU2), we also found the 
minus contribution, uxd^ui — u^d^ux. The plus contribution is sourced by both the first 
order off-diagonals and by the second order diag onals. The minus term plays an important 
role in the approach advocated in Refs. Ill ij especially near the degeneracy (small A), 
where it exhibits a resonant enhancement. When compared with our results, in the region of 
near degeneracy we find a weak enhancement in all contributions to the CP-violating vector 
current, such that the minus contribution remains subdominant. 

By performing a numerical study of fluid equations, we have analysed the CP-violating 
vector and axial vector currents for a slowly moving wall, and found that, in the nonlocal 
regime, in which the currents are only weakly damped TL W < 1, the first order terms provide 
a dominant contribution to the CP- violating currents if A = mf—m 2 (m? (i = 1, .., N) denote 
the mass eigenvalues) is smaller than about 20 T 2 . Whether this statement remains true with 
respect to the baryon asymmetry remains unclear, since the oscillations, poor transport and 
the tracelessness of the first order terms could prevent an efficient production of BAU. 

A more comprehensive comparison with the work of Refs. Q] and jloL IrI \v\ . and 
the explicit calculation of the BAU is postponed to a forthcoming publication. 

Acknowledgments 

We would like to thank Steffen Weinstock and Marco Seco for numerous useful discussions. 

APPENDIX A: DIAGONALIZATION OF THE CHARGINO-HIGGSINO MASS 

MATRIX 

The chargino-higgsino mass matrix is given by 

rnJ M ^ Hi ) (Al) 

\gH"i ^ } 

The mass matrix m is diagonalized by the biunitary transformation 

m d = UmV ] , (A2) 
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with 

V=(^Y( UA + A) a I (A3) 
U(A + A); ^ _ a * I( A + A) ' 

a = g(M 2 H 1 + y* c H*) , A = |M 2 | 2 - |/i c | 2 - (u\ - u\) , A = (A 2 + 4|a| 2 )^ 

and 

VAIA + A); ^ |(A + A) J 1 * 

a = g(M*H* + y^Hx) , A = |M 2 | 2 - |/i c | 2 + {u\ - u\) , A = (A 2 + 4|a| 2 )5 = A 

where we defined U\p = \gHxp\- Note that a and A can be obtained from a and A by the 
replacements, M 2 <-> M 2 , // c A** an d i?i <-> H\, such that A = A, as indicated in (|A4|) . 
The mass eigenvalues-squared are given by 

™*i/2 = \{\M2\ 2 + |/i c | 2 + (u? + u\)) ± | (A5) 
and can be calculated quite simply by noting that UmrrvW = m? d = VrrvmV* . 
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